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Abstract 

We examine the growth of crystals from vapor. We assume that the Wulff shape is 
a prism with a hexagonal base. The Gibbs-Thomson correction on the crystal surface is 
included in the model. Assuming that the-initial crystal has an admissible shape we show 
local in time existence of solutions. 
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1 Introduction 

The main result of this paper is a mathematical study of crystals grown from supersaturated 
vapor. We assume that evolving crystal Q is a prism with A-gonal base. We assume that 
the number N is constant. Let us stress that our crystal fl has not to be a convex set (see 
Figure 1 in the case N = 54). Such kind of ice crystals are formed in the atmosphere. Our 
activity is motivated by Gonda and Gomi results, [H], and also by Giga and Rybka [9]. 
The authors of the paper [9j assumed that the Wulff shape is a fixed cylinder and also that 
the process was slow, i.e. they considered quasi-steady approximation of the diffusion equation 
(diffusion is much faster than the evolution of free boundary) . Giga and Rybka obtained results 
for quasi-steady approximation system (see [TO], [UJ, [32], [13]). 




In this paper we do not assume that the crystal evolves slowly. Namely, we assume that su- 
persaturation a fulfills diffusion equation with a drift 

e^a = Aa + f-Va (1) 

outside of crystal Q. For simplicity of notation we assume that e — 1. We are not going 
to examine the properties of the solution when we change the parameter e. In the above 
equation / is a velocity of vapor. We assume that / is given. In fact, velocity / should satisfy 
an addition equation. 

We assume that the supersaturation a has a specific value at infinity, i.e. 



lim a(x) = a°°. (2) 

♦oo 



We require that mass conservation law on the crystal surface is fulfilled. Namely, 

da , , 

where v is the outer normal. This condition is very natural in the epitaxy model. See Eliot and 
Gurtin paper [6j. Let us stress that it could be better to discussed the system with replaced 
g by V in ([3D, where V is the velocity of the growing crystal. Namely, V = This is the so- 
called Stefan condition. Our method can not be applied to this case. This problem will be 
examined in the forthcoming article |15| . 

The value of the supersaturation a at the surface satisfies Gibbs-Thomson law, namely, 

- a = -div s £ - (3V, (4) 

where (3 is the kinetic coefficient and £ is a Cahn-Hoffman vector field (see [IE])- We denote by 
V the velocity of the growing crystal. We also assume that the velocity V is constant on each 
side of the crystal. Our model does not include bending and breaking of surface dQ. Roughly 
speaking the above relation says the supersaturation on the crystal surface is proportional 
to the curvature of surface (— divs£)i an d to the velocity of the evolving crystal. Let us discuss 
the properties of the Cahn-Hoffman vector field £. Let 7 be surface energy density (Finsler 
metric, see [27]), i.e. 7 : M 3 — > M is a 1-homogeneous, convex and Lipschitz continuous function. 
If the surface of the crystal dfl and the surface energy density 7 are smooth, then 

= V 7 Mx)). (5) 

If we do not assume that 7 and dfl are smooth then the above expression may not make sense 
(see [8]). Instead V7 we can write subdifferential dj, which is defined everywhere. In this case 
c>7 is a nonempty convex set (not necessarily singleton). The condition ((SJ) can be replaced by 
the following relation 

£(x) G d r ){p{x)). 

We are not going to discuss the properties of the Cahn-Hoffman vector field. In preliminaries 
section we show a "nice" property of the averaged divergence of the Cahn-Hoffman field £, 
namely 

—KiH 2 {Si)= I divsid7i 2 {x) ) 
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where is a crystalline curvature (see next section). Let us stress that the LHS of the above 
expression is independent of £. 

In this paper we assume that the Frank diagram F 1 (see pEE] or next section) is a sum 
of two regular pyramids having a common base. Hence, the Wulff shape W-y is a prism with 
a hexagonal base. We assume that the crystal is an admissible shape, i.e. the set of outer 
normal vectors to dfl coincides with the set of outer normal vectors to W 7 . Let us notice that 
Q need not be a convex set. Let us stress that we could apply our analysis to more general 
Wulff shape. Namely, we could assume that W 7 is a prism with regular fc-gonal base, where 
k > 3. But only hexagonal symmetry is intersting from pysical point of view (see Figure 1). 

In order to simplify our system (Q]), (El), ([!]), (J4]) we will consider averaged Gibbs-Thomson 
law. Taking into account the above considerations and the assumption that the velocity V 
is constant on each side of the crystal we obtain 



Let us mention that averaged Gibbs-Thomson is not new condition. Namely, this equation has 
appeared in [19] and [2l] . 

Our main aim is to show the local in time existence of solutions to (pQ), ([2]), Q, ([6]). We show 
this in a few steps. Let us notice that we can look at this system as parabolic equation coupled 
with ordinary differential equation for the evolving sides of crystal (the signed distance side 
from Si(t) to Si(0)). Next, we transform our problem to the system in the fixed set, in this case 
it is outside the initial crystal. In this domain we solve our system. Subsequently, we examine 
the linear problem and here we apply the analytic semigroup theory (see [I], [2], [4], [7], [20] . 
[23]). Next we construct an approximating sequence and apply Ascoli's Theorem. Let us notice 
that we work in an unbounded domain and the standard Rellich-Kondrachov compactness 
theorem does not work. We have to show additional properties of the approximating sequence. 
Our method is similar to the method presented in [28] . 

The paper is organized as follows. First, we introduce notations, recall definitions, show 
some facts about crystalline curvature and formulate the problem. Next, we show the main 
result. The last section (Appendix) contains useful result. 

At the end of this section, let us comment on the mathematical literature. In the smooth 
case, i.e. assuming that the surface of the crystal is differentiable manifold and for the two 
phase Stefan problem this model was solved in the class of smooth function (see Chen and 
Reitich [5]). Independently, Radkevich (see [26]) showed local in time existence of a smooth 
solution to the problem considered by Chen and Reitich. The advantage of Radkevich is that 
the author allows a slightly more general form of the diffusion equation. The problem for ft = 
and smooth interfaces was studied by Luckhaus [22] and in greater generality by Almgren- 
Wang [3J. In particular, they showed that uniqueness fails. Let us stress that these authors 
have worked in bounded domain. Quasi-steady approximation in the case when the crystal is 
a cylinder and a domain is unbounded was discussed by Giga and Rybka [9J. Let us mention 
that these authors discussed the system with replaced g by V in (jHJ), i.e. the so-called Stefan 
condition. Our method can not be applied to this case. 

Notation. We use the following convention, whenever we see the inequality A < cB we 
tacitly understand that it holds with some positive constant c independent of A and B. 




(6) 
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2 Preliminaries 



Let us denote an evolving crystal by fl(t), and its exterior by Q(t) c = M, 3 \Q(t). Let S(t) = dQ(t) 
its surface. We assume that crystal Q(t) is a prism at all times. To be more precise we write 

Sl(t) = {(x 1: x 2 ,x 3 ) G R 3 : L B (t) < x 3 < L T (t), (x 1: x 2 ) G W(t)}. 

Where W(t) is an iV-gon in the plane. 

Let us divide the surface into N + 2 pieces, i.e. 

Si(t) = {(xi,x 2 ,x 3 ) : (xi,x 2 ) e Li(t),L B (t) < x 3 < L T (t)} , 
S T (t) = {(x u x 2 ,x 3 ) : (x u x 2 ) e W(t),x 3 = L T (t)}, 
Ssit) = {(x!,x 2 ,x 3 ) : (x!,x 2 ) e W(t),x 3 = L B (t)}, 

where % = 1, 2, N (we also use the notation N + 1 = T and N + 2 = B). The sets Si shall 
be called facets. In order to understand this convention it is better to look at Figure 2. 



Figure 2 

Let us denote by L(t) the boundary in the plane of 7V-gon W(t), namely 

N 

L(t) = dW(t) = (jLi(t), 

i=l 

where Li(t) are edges. Let us additionally write Qq = Q(0) , and the coordinate systems will 
be presented as follows 

(t,x 1 ,x 2 ,x 3 ) G R + x Q c 
(t,xi,f 2 ,x 3 ) G {t} x tt c (t). 

We denote by n the outer normal vector to the boundary Q$ and by v the outer normal 
to tt c (t). 

We denote the coordinates of vertices of Qo by 
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where % = 1, N and L is the height of Qq. Through the paper we shall assume that G O(0). 
Analogously, we denote the coordinates of the vertices of O(t), i.e. 

(x\(t),x\(t),L T (t)) , (x{,xi,L B (t)) , where i = l,...,N. 

Next, we denote by Zi(t) the signed distance between i-th face of O(0) and z-th face of O(t), 
i.e. 

(dist{Si(t),Si(0)), if (xi(t)-4(0),x J 2 (t)-xU0),L T (t)-L T (0))-z/ 4 > 0, 
" " 1-^(^(^,5,(0)), if (xi(t)-xl(0),x J 2 (t)-x^(0),L T (t)-L T (0))-z/ l <0. 

Making not difficult but boring and long calculations one can examine what is the dependence 
between Zi(t), Xj(t) and X/~. We will summarize them. We use the following convention Xj = 
Xj(t). Namely, let us denote by 

Cj — ^2 ^2 5 ^i — X~\ (0) 

di = Xi x'i , fi = x 2 (0) x 2 . 



Next, we denote by 



fcj — 

We can assume that x\ ^ for i = 1, N. Hence, the normal vector has the form 



ru 



\.L 1 X 2 | 



-,1 

di 



Via) 2 + 

Finally, we can state the following lemma. 
Lemma 1. 

ll ' ~ m (t) ' 21 ' ~ m (t) ' 

where 

xiai(t) = (x\ + ai_iZj_i)(cj_i + 6j_i^_i)(dj + a,iZi) - (x\ +1 + aj^)(cj + b i z i )(d i - 1 + aj_i^_i) + 
-(x 2 + 6 i _i2! i _i)(d i + aj2;i)(dj_i + a^Zi-i) + (x 2 +1 + 6»^i) (d» + a^X^-i + aj-i^-i), 
a?2ai(*) = 0*a + ai_i^_i)(cj + 6^) (<;»_! + 6 i _ 1 2 i _ 1 ) - (4+ 1 + a i ^ i )(c i + b,^) (ci_i + + 
+ (x 2 +1 + biZi)(di + aiZi)(ci-i + fc^i^-i) - (x 2 + ^-i^-i) (c?i_! + a i _i2 i _ 1 )(c i + 6^), 

= (ci-i + + a,Zi) - (q + biZi)(di-i + di-iZi-i). 

Proof. We leave those straightforward calculations to the interested reader. 

□ 

In our problem evolution is determined by normal vectors to the boundary of the initial 
crystal. It means that facets can move parallel to the facets of O . We assume that the set 
of normal vectors of O coincides with the set of normal vectors of a prism with a hexagonal 
base. Such O we shall call an admissible shape. Let us notice that we need N + 2 parameters, 
in order to describe the evolution of crystals, i.e. distance between facets Z{. In our problem 
there appears surface divergence divsC- We define this quantity as follows 

divs£,(x) := tr((Id — n(x) ® n(x))V£(x)), 

where n(x) is the unit vector orthogonal to tangent space T X S (see |29|). 
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2.1 Crystalline curvature and surface energy 

Let us recall the basic objects from crystalline geometry, i.e. Wulff shape W y and Frank 
diagram F 7 . Let 7 be a surface energy density, i.e. 7 : R 3 — ► R is a 1-homogeneous, convex and 
Lipschitz continuous function. Then the Wulff shape and Frank diagram (see [18] ) are defined 
as follows: 



In our paper we assume that the Wulff crystal W 1 is a prism with a hexagonal base. This as- 
sumption is consistent with physical experiments (see [25] and Figure 1). Then, it is not hard 
to see that the Frank diagram F 1 is a sum of two regular pyramids having a common base. 
Let us formulate the following 

Lemma 2. Let us assume T 6 G 5*0(2), i.e. 



where Ai is i-th sector (i = 1, 2, ...Q), see Figure 3. Then F 1 is a sum of two regular pyramid. 



W-y = {1 £ R 3 : Vn G M 3 , \n\ = l,x ■ n < 7(71)} 
F 7 = {x G R 3 : 7(z) < 1}. 




and 7(1, 0, 0) — Jl- Let us also assume T 6 (l, 0) T = (x' 1; x' 2 ) T and 




and 



7(21, 0:2,23) = 7i (^6 (i 15 (^1.^2) ) + \x3\lT if (x!,x 2 ) G A 




Figure 3 
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Proof. In this relatively easy calculation homogenity of function 7 and elementary geometry 
have been used. The details we left to the reader. □ 

Let us recall that the surface energy E(S) is expressed by the formula: 

E(S) = [ j{n{x))dH 2 {x), 
Js 

where H 2 is two dimensional Hausdorff measure. 

The crystalline curvature «j is defined as (see |32j): 

Ki = - lim — ^, 

a^O AV 

where a is the amount of motion of Si in the direction of the outer normal to S, AE is the re- 
sulting change of surface energy, and AV is the change of volume. Let us denote by 014 the angle 
between z-th and % + 1-th facet. Now, we show the following lemma. 

Lemma 3. The crystalline curvature is given by expressions: 

«* = -faM + ctg{a J+l ))l{n,) - ^tK-i) - ^-^7^+0 + yTK) 
Where lj = \Lj\ and I = \L\, j = 1,...,N. 

Proof. It is a straightforward calculation. It is not hard to see that: 

N 

AE T = ^Ija-finj) 
i=i 

AV T = aH 2 (S T ) 
AEj = la\cbg{aij) + ctg(aj + i)jj(nj)+ 

+ ■ 1<X ( + ■ / " , ,M n j+1 ) + a(2lj + a{ctg{a j ) + ctg(a j+1 )))^(n T ) 

sin{ot.j) stn{ctj + 1) V v / 

AVj = i( 2 ^' + a ( ct9 ( a ^ + ct 9{ a j+i))) i i 

where j = 1, N. From those formulas the claim follows. □ 

We finish this section with a theorem which let us average Gibbs-Thompson law. 
Theorem 1. The following equality is true: 

1 



divs^dH 2 (x). 



Proof. It is a similar calculation to the proof of Proposition 2.1 from [12], where the Gauss 
formula was the main tool. We leave this to the reader. □ 
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Finally, we will work with the system 



^-a = Aa + / ■ Va in fi c (t) 
dt J w 

g = ^ on dfi c (t) (7) 
o{x)H 2 {x) = -Ki^+PiVi, 



where g = (g x , ...,g N+2 ) is given. 

3 The existence of a solution - preparation 

Let us mention that we can assume that a°° = 0. At the beginning we transform the system 
to the one with homogeneous boundary conditions. From Lemma fTUl (see Appendix) we know 
that there exists a family of functions hi, where % = 1,2..., N + 2 with compact supports, such 
that 



dhj 



= Si j. 

5,(0) 



Subsequently, let us introduce the notation 

JV+2 



9i 



i=i 



and u = a — G. Now we obtain the system for u with homogeneous boundary conditions. 
Namely, we will deal with the following final problem 

— = Au + f -Vu+ (A + f -V - - ) G in Q c (t), 



fit i 

^ = on dQ c (t), (8) 
av 



1 



n - 1{ , , {f)) ^ t (u(x) + G) dH 2 (x) = -Ki(t) + m, 

u(0) = u , z{0) = 0. 



We assume that / G L c 

3.1 Transformation of the system to the fixed domain 

In this subsection we transform our problem into a system in a fixed domain. In this case it 
is the outside of the initial crystal Q . We have to construct a family of diffeomorphisms from 
n c to tt c (t), i.e. 

A : Q c -> Q c (t). 

Let us notice, that Q c (t) depends on z(t), namely Q c (t) = Q c (z(t)) = Q c (z). We have the fol- 
lowing main theorem of this section. 
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Theorem 2. There exists a family of diffeomorphisms A z of C 2+a class (a > 0) trans- 
forming Q c (0) on Q c (z). Additionally, the family has the following property, if \xi,x 2 ,x 3 \ > 
10 diam(Q c (z)) then A z = Id. In addition, the dependence upon z is smooth. 

Before we go to the proof we make some comments. 

Remark 1. // we make a cross-section parallel to St, then one of the problems is to find 
a diffeomorphism transforming the exterior of N-gon onto exterior another N-gon. It seems 
that we could use complex analysis methods and apply Riemann Theorem. But here we meet 
the following problem, namely we work in nonsimply connected domain. There exist the so-called 
Christoffel-Schwarz expressions transforming half plane onto interior of N-gon. The second 
problem is the requirement of smoothness of A with respect to t. 

Remark 2. The number 10 in formulation of the theorem is not substantial. We could take 
any number greater than 1. 

Remark 3. During the examination of the transformed system we use the property that A z is 
identity far from S(z). 

Proof. We divide the construction into two steps, namely 

$ : Q c -> E(z) 
V : £(*) -» n c (z), 



where £(z) is defined as follows: 

£(*) = {(x u x 2 ,x 3 ) G M 3 : L B (z) <x 3 < L T (z) 
We define the map A in the following way 

A = \fr o $. 



;r 



i,x 2 )eW(0)} 1 



First we define transformation $ as follows 



(x 1 ,x 2 ,L T (z) + (x 3 - f ) f 

(xi,x 2 ,x 3 + L T (z)) 
(xi,x 2 ,x 3 - L B {z)) 
(x 1 ,x 2 ,(f) 1 (x 3 ) 
(x 1 ,x 2 ,4> 2 (x 3 ) 



L T (z)-L B (z) 
L 



if x 3 E 
if x 3 G 

if x 3 e 
if x 3 e 

if x 3 G 



f + 25, 00) 
-00, -| -25] 
± + 5,%+25] 
-±-25,-±-5] 



where 5 > is a fixed number and 4>i,4> 2 are smooth and one-to-one functions such that 
$ is smooth. Let us notice that $ for x 3 close to facets is scaling in direction x 3 . 

Now, we start to construct the map The idea comes from differential geometry and is sim- 
ilar to exponential map exp (see also construction of Lie algebra q of Lie group G), we refer 
interested reader to the textbook [3 1 J . We construct a vector field T which transforms the 
domain £(2) onto fl c (z). Next we integrate this vector field and we obtain the flow. This flow 
gives us our diffeomorphism. 

Let us define the smooth map f] 



n 



[0,1]. 
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Figure 4: 



We require that 



r)(x 1 ,x 2 ) 



1 if (x!,x 2 ) £ W(0) or dist((x 1 ,x 2 ),dW(0)) < § 
if lx!,x 2 ) e W(0) and dist((x 1 ,x 2 ),dW(0)) > 



2s 

3 ' 



We will need it below. 

Subsequently, let us define the map t as follows 



t(xi,x 2 ) 



a' i (x 1 ,x 2 ) 



a' i (x 1 ,x 2 ) + a i {x l ,x 2 y 
where ctj and a[ are explained in Figure 4. 



Next let us take a smooth map: 

X- [o,i] -[0,1], 

such that X[o,|) = 0, = 1 an d x'(0) — x'(l) = 0- Let us define r in the following way 

t = x°t. Let us notice that we can not take simply r = t. Our vector field would have required 
the properties without lines determined by directions Wi. 
Subsequently, we define the map & as follows 

& = (1 - r(xi,a;2))wj + t(x u x 2 )w i+1 . 

Our diffeomorphism must fulfil the condition: the boundary of iV-gon is transformed onto 
boundary of iV-gon. Now, we have to construct the vector filed T in a smooth way. Namely, 

{6 if (xi,x 2 ) e Ai(z) and x 3 e (L B (z), L T (z)), 

■n(xi,x 2 )€i if (x 1: x 2 ) E Ai(z) UOi(e) and x 3 ^ (L B (t), L T (z)), 
in another cases. 

Where Aj(z) is i-th sector (see Figure 5) and 

O t (e) = (L t (0)) £ \(a t (z)Ub l (z)) 

where (A) £ is the e - envelope of set A, i.e. 

(A) £ = {x : dist(x, A) < e}. 
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Figure 5: 



Si(z) 




3i(z) 



bi(z) 



Figure 6: 



The sets a,i(z) and bi(z) are explained in Figure 6. Thus we have obtained a smooth vector 
field. Integrating this vector field we obtain the flow T t (xi,x 2 ,x 3 ). We define transformation 
\& as follows 



From the proof of the above theorem we obtain the following 

Corollary 1. Let J be the absolute value of Jacobian of transformation A, then: 

a) J e L^, 

b) if \zi\ < M , where i — 1, N + 2, then there exist /, L such that < I < J < L, 

c) A is a function of (xi, x 2 , X3) and Zi(t). 




Finally our diffeomorphism has the form 




where $3 is the third component of $. 



□ 
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Now, we write the system (jSj) in the new variables. It is not hard to see that 



(xi,/r 2 ,x 3 ) 



d_ d_ 



{^(x 1 ,x 2 ,x 3 )A(x 1 ,X 2 ,X 3 ,t)) V 

(OA \ ^ 
t^t) (V( Xl , X2 , X3 )A(x 1 ,x 2 ,x 3 ,T)y 1 X7 



(xi,X2,X3) 



(V( xuX2>X3 )A(x 1 ,x 2 ,x 3 ,t)) 1 n 
We will use the following notation: 

u(xi, x 2 , x 3 , r) := u(x~i, x 2 , x 3 , t)\^ xliX2jX3> t^^( xljX2 ,x3,T),r)- 
Let us denote by H(t) the transformed Laplace operator and 

H(t,z,z)V( XuX2jX3 ) = \ +f (V ) A(x 1 ,x 2 ,x 3 ,t)) V 

The boundary operator is defined as follows 



5 (t)a 



(V (a;iiX2iX 3)A(xi,X2,X3,r)) *n 
where A, B and G are defined in the following way 



(A(t),B(t),C(t)) 



a 



Mr) = ^ (A 2 x Al 3 d Xl - AlA 3 X3 dx 2 ) 

B{t) = i - 

C(t) = ^ ((Ai 2 A^ - «) X1 + (A^ - A^Ay ^ + (A^ - A^J d X3 ) 

Q ~ A X3 (A Xl A X2 — A X2 A X1/ 

We have used the notation A l x . = Finally, system (jSJ) takes the form 

= H{z(t))u + ff(-z(t), i(t)) Vu + ^ff(t) + #(t)V - j G in fig, 

B (t)u = on 8fig, 

— / ( M (x) + G)j(x, 2 ;(t))^ 2 (x) = «: i (t)+AV r i , 

u(0) = Uq, -2(0) = 0. 

where /, G, u are old /, G, u in new variables. 
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4 The existence of a solution - a linear theory 

In this section we show solvability of the following problem 

d 

—u = H(t)u -5u + f(t) in fig, 

B (t)u = on dfig, 
u(0) = M , 

where 5 will be specified later and / is a Holder continuous map. 
Let us define an operator with the domain 

H(t) : D(H(t),Q c ) -> L 2 (Q C ) 
D(H(t),n c ) := {v G F^fi ) : #(t)v G L 2 (O c ), S (*)«|fln«= = 0} . 

Where the norm is defined as follows 

\\u\\D(H(t),n°) ■= (\\u\\%i(nc) + \\H(t)u\\ 2 L 2^ c) ^ . 
Later we will need the following 

Proposition 1. D(H(t),Q c ) is Banach space for each t > 0. 

Proof. We know that D(H(t),Q c ) and D(A,Q c (t)) are isomorphic (see the proof of Corollary 
02). Hence, it is enough to show that _D(A,fi c (i)) is Banach space. 

Let us take Cauchy sequence {uk} from D(A,Q c (t)). Subsequently, D x .Uk and Aw fc are 
Cauchy sequences in L 2 (Q c (t)). So D x .uu — > v i: and Aw fc — > / in L 2 (fi c (t)). In particular 
ma,. — ► it in L 2 (VL c (t)). Hence, 



uD Xi <pdx <— / UkD Xi (f)dx = (-1) / D Xi Uk4>dx — ► (-1) / t^dar. 
n c (t) ifi c (t) in c (t) Jf7 c (t) 



Analogously 



uAcfidx <— / UkAcpdx = (—1) / Auk4>dx — ► (—1) / fcj)dx. 
n c (t) Jn c (t) Jw(t) Jn c {t) 

Hence, we obtain -D^w = t>j and Am = /. From trace theorem and D(A, fi c (t)) C if 2+ £ (fi c (t)) 
(see [16]) we obtain 



<9(w; - it*) 



< C||lij — Mfc||z)(A,Q c (i))- 

L 2 (9Q c (t)) 



<9z/ 

Hence, we obtain the assertion of the proposition. □ 

Now, we show the following 
Lemma 4. If u G D(A,fi) then there exists a sequence {R n }, Rn — > oo such that 

lim / -^u dU 2 = 
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Proof. From the integration by parts formula the following expression follows 
/ uAudx + / \Vu\ 2 dx = ( —udTi 2 . 

Jn c \B(0,R) Jn c \B(0,R) JdB(0,R) d v 

But it, Vh, Au G L 2 (Q c ). Hence there exists a sequence R n such that 

uAudx — > 0, 

Q c \B(0,R n ) 



\Vu\ 2 dx — > 

Q c \B(0,R n ) 

as R n — >■ oo. 

From this the lemma follows. □ 

Now, we can go to the following result. 

Theorem 3. Each of the operators —H(t), t>0, generates an analytic semigroup. 

Proof. It is enough to show that —H{t) is a sectorial operator (see [23] or [20]). First of all 
we examine the operator H(0), i.e. A. Namely, we show the following, 

Lemma 5. The operator —A with the domain D(A,Q C ) is sectorial in L 2 (fl c ). 

Proof. First of all we show that the operator —A is closed and densely defined. Indeed, 
since C °°(^) C £>(A,fi) and C °°(fi) is dense in L 2 {Vt), hence D(A,Q) is dense in L 2 (fi). 
Now, we show that —A is a closed operator. 
We need the following estimate 

l|Vv||i»(nc) < IMlL2(f^)||Aw|| L 2 (nc) (10) 

for v G D(A,Q C ). In order to show inequality f fTOl) we take the following equality 

Av = f 

where / G L 2 (fl c ) and v G -D(A). We multiply both sides of the above equality by v and after 
integration by parts we get, 

ll^lli 2 ^) = / v fdx < ||w|U 2 (Q c )||/IU 2 (n c ) = IMU 2 (Q c )||At;|| L 2 (QC) , 

where we used Lemma [4] and Schwarz inequality. 

Now, let us take a sequence {t^}^ C D(A,Q C ) such that: 

v k -> v 
Av k ^f 
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in L 2 (Q C ). From inequality ( fTOl) we obtain 



\\ v k ~ Vi £)(A,f2 c ) — \\Vk 


- Vi\ 


|2 


+ ||Au fc 


- Au,| 


|2 

lL 2 (fi c ) — 


\\vk ~ vi\\h(n°) + llV^fc - 




|2 

lL 2 (n c ) 


+ ||Au fc 


-Avi\ 


|2 < 


IK ~ v i\\h(nc) + \\ v k ~ vi\\ L 2(nc)\\Av k - 


Au,| 


|i 2 (n c ) 


+ ||Au fc 


-Avi\ 


|2 < 


(Ik 




L 2 (n c ) 


+ ||Au fc 


-Avi\ 


Il 2 (^)) 2 



Hence 

Ikfe - ^lb(A,n=) < \\vk - ^IU 2 (n e ) + 1 1 Aw* - A^|| Z 2 (nc ). 

Subsequently, we obtain that {vkj^Li is Cauchy sequence in D(A,Q C ). Finally, 

v k -»• v w L>(A,fi c ). 

Hence, t> G _D(A,f2 c ) and Af fc — > Aw. 

Now, we show the estimate for a resolvent operator. We follow the method used in the 
book [20]. Let us consider the equation 

Xu + Au = f 

where / G L 2 and u G D(A). Let us assume that ReX < 0. 
Let us multiply the above equation by u and integrate: 

A|w| 2 + uAu = I uf. 

Integrating by parts we obtain 

' AM 2 - |Vu| 2 = / uf, 



where we applied Lemma [H 

Next, we take the real and imaginary part of above expression: 



i2eA||tt||| 2 — ||Vm|| 2 2 = Re uf 

Jn c 

— Im\\\u\\ 2 L 2 = Im / uf. 



Hence, 



-ReA||M|| L 2 (nc) < ||/||i2 (n< 



\Im\\ ||U|| L 2( QC ) < \\J ||L2(Qc). 

Subsequently, we take a sum of the above inequalities 

|A| ||w|U 2 (n c ) < 2||/||l2(qc). 
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Next, for arbitrary A we behave in a standard way. It means, if ReX > and < | arg(A)| < tt, 
let us denote A = a + ib we obtain tan0 < |^|. From the above estimates we obtain: 

|6| ll M IU 2 (nc) < Il/||x2(n<=) 

and 



^0 



Hence, 



|A| |M|£2(nc) < ^^ll/lli 2 (^)- 



This ends the proof of the lemma. □ 
Now, we show Theorem [3] for arbitrary t > 0. 

Using the same argument as proof of Lemma [5] one can show —H{t) is densely defined. 
Subsequently, we show that the operator —H(t) is closed. First of all we show the following 
inequalities 

\\^ u \\h(nc) < c\\H(t)u\\ L 2 {nc) \\u\\ L 2 {nc) , (11) 
where u G D(H(t),Q c ). Indeed, let us consider the equation 

H{t)u = f in Vl c . 

Changing the variables in the above equation we obtain the following one 

Am = / in n c (t). 

where u and / are old u and / in the new variables. Following the method presented in 
the proof of Lemma [5] and from Corollary 1 we obtain the inequality ( fTTl) . 

From inequality (fTTTl closedness easily follows. Indeed, let us take any sequence {uk} C 
D(H(t),Q) such that Uk — >• u and H{t)uk — ► / in L 2 (fl). Hence Uk and H(t)uk are Cauchy 
sequences in L 2 (fl). Subsequently, from inequality (PTT1) we obtain the following estimate 

||V«fc - V« z ||| 2 ( n c) < c||if(t)(n fc - w/)|| L 2(nc) ||wfe - «i||x2(n<=) 

and this inequality finishes the closedness of the operator H(t) (see proof the Lemma [5] for 
details) . 

Now, we show the estimate for the resolvent operator. Let us consider the equation 

Am + H{t)u = f in fi c . 

where / G L 2 and w G D(H(t),Q c ). Changing the variables in the above equation we obtain 
the following equation 

Xu + Au = f in n c (t), 
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where u and / are old u and / in the new variables. 

From the proof of the Lemma [5] we obtain the inequality 

|A| |M| L 2 (n c) < C||/||i2( n o). 

After changing the variables we obtain 

|A| I \u\ 2 Jdx <C I \f\ 2 Jdx. 

From Corollary 1 the following 

\M IMU 2 (n c ) < Cj\\f\\ L 2 (n c) 

holds. 

From this the theorem follows. □ 

Before we state next result we recall the definition of some hypothesis (see p] and |2j). 
Namely, 

Hypothesis I. For each t G [0, T], A(t) : D(A(t)) C I? — > I? is a closed linear operator 
and there exist M > and 9q G (|, ir) such that 

p(A(t)) D S eo = {z 6 C : | arg^| < 6 } Vt G [0,T], 
||i?(A, A(t))|| £(L2) < VA G S eo ,W G [0,21. 

Hypothesis II. There exist B > 0, k G N, ai, ...,atk, Pi, •••■> Pk with < /3j < a« < 2, such 
that 

k 

1 1 J2(A, [Ait)- 1 - A(s)- 1 ] \\c(l 2 ) <B^2(t- sr 

i=l 

VA G S 6q \ {0}, V0 < s < t < T. Where R(X, A(t)) is resolvent operator. 
We have the following 

Theorem 4. There exists 5 > suc/i iftai i/ie operator 

A(t) = H(t) - Sid 
with a domain D(A(t)) = D(H(t)) fulfills Hypothesis I and Hypothesis II. 

Proof. The proof of Hypothesis I is almost the same as the proof of Theorem [3l 
Concerning Hypothesis II, if we set for/ G L 2 (Q C ) 

v = A(s)- 1 /, u = R(X, A{t)) [X - A{s)]v, 

then we have to estimate L 2 -norm of 

it — v — A(t)R(X, A{t)) [Ait)' 1 - A(s)- 1 ] f. 
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Now, u and v solve respectively, 

Xu - A(t)u = Xv - / in Vl c 
B(t)u = on dVl c 

and 

A(s)v = -f in n c 
B(s)v = on dVl c . 

Hence, u — v solves 

X{u -v)+ A(s)v - A(t)u = in Q c 
B(t)u = B(s)v on dVl c . 

Next, using the similar methods as in the proof of Theorem [3] and from inequality 

||V«|| £ 2 (n c) < c||/|| X 2 (n c) 

and Lipschitz continuity of coefficients of the operator A the proof follows. □ 

Now, from the above theorem and [I] the following holds. 

Corollary 2. If f is a Holder continuous map, then there exists unique solution of the following 
problem 

d 

—u = A(t)u + f(t) onn c , 
B (t)u = in dn c , 

U(0) = M - 

It is given by the formula 

u(t) = U(t,0)u + [ U(t,a)f(a)da. 
Jo 

5 Nonlinear theory 

Before we state the main result we are going to present some facts. 

In order to understand better interpolation spaces we show isomorphism theorem. 

Corollary 3. The following interpolation spaces are isomorphic 

[L a (n& D(A, nt)] e - [L 2 (Q C ), D(H(t), Q c )] e . 
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Proof. From Theorem [2] we know that there is a family of C 2+a diffeomorphisms A t : Vtl 
Let us denote by A t „ the pullback operator. Hence we can write the following diagram 



L 2 m 




At, 



(L 2 (^), D(A, ^ (L 2 (fi c ), D(H(t), n c )) 



■ L 2 (n c ) 




D(H(t),tt c 



A, 



So we obtain an isomorphism on the level of interpolation pairs. It is well known (see |33| ) 
that interpolation is a functor from a category of Banach pairs into a category of interpolation 
spaces. Hence, from functoriality (see [30] or [24]) we obtain 



(L 2 (fi t c ),£>(A,fi t c )) 



[L 2 (n?),D(A ) n?))] t 

Where 6 is the interpolation functor. 



(L\n c ),D(H(t),n c )) 



[(L 2 (fi c ), D(H(t), Q c ))]( 



□ 



Corollary 4. // (3 > 



e ?5 sufficiently small and s < | + e, i/zen t/ie following is true 
[L 2 (Q c ),D(H(t),Q c )] ^ H S (Q C ). 



Proof. It follows immediately from the properties of interpolation spaces and Grisvard theory 
(see [Ej and [17]), see also [33]. □ 

Now, we show a compactness result. 

Lemma 6. Let us assume Q C IR fc is bounded domain with Lipschitz boundary. Let us also 
assume that the sequence u n is bounded in H cr (R k \ Q) and | |itn| |.ff CT (R fe \,B(o,7V)) ~ * as N oo 
uniformly with respect to n. Then there exists subsequence of u n convergence in H a (M. k \ Q) 
for a < a. 

Proof. We can assume that Q C 5(0,1). Let {u n } C H a {R k \ Q), then {u n \ B (o,N)} C 
H a (B(0, N)\Q). Using Rellich - Kondrachov Theorem we can extract a sequence {w n i|s(o,i)} c 
H a (B(0, 1) fl fl) such that {lim I b(o,i)} ~~ ¥ u i m H a (B(0, 1) \ O). Let us denote by u\ exten- 
sion of U\ onto M fc \ fl. So we have chosen a subsequence {u ni } (we restrict the sequence 
to such that the following inequality is satisfied \\u m — &i||.ff<*(.B(o,i)\fi) < 1)- Analogously 
we work on ball 5(0,2), e.g. we subtract subsequence {w n i|.B(o,2)}, (we will denote it by 



19 



{ u n 2 1 .8(0,2)}) such that {u n2 \B(o,2)} u 2 . Denote extension by u 2 and we take only this terms 

such ||w n2 - u 2 \\h"(b(o,2)\q) < \ end so on, namely 

u ni — * Ui in H a (B(0, 1) \ Q) and extension U\ 

and such that ||w ni - ui\\H<*(B(o,i)\n) < 1 

u n , 2 — > u 2 in H a (B(0, 2) \ f2) and extension u 2 

and such that ||w„ 2 - U2|if«(s(o,2)\n) < | 

u n3 — > M3 in H a (B(0, 3) \ f2) and extension -u 3 

and such that \\u m - W3||i?«(B(o,3)\Q) < § 



Mnjv mat in H a (B(0, N) \ Q) and extension mat 
and such that \\u nN - ujv||H«(8(o,JV)\n) < ^ 

where w njv is subsequence of «„ M for iV > M. We shall show that {un} is convergent in 
H a (R k \ Q). Indeed, we show that Cauchy condition is fulfilled 

\\UN — UM\\H a (R k \n) < \\Un — UM\\H<*(tlnB(0,min{N,M}) + 
\\^ L N\\H a («. k \B(0,min{N,M})) + \\^M\\H a (R k \B(0,min{N,M}))- 

Hence, we obtain 

limsup||MAr - u M \\H"(Rk\n) = 0. 

N,M^oo 

Let us denote the limit of this sequence by u. We show that u nN converges to u in H a (M, k \Q). 
Namely, 

\\Un N — u\\ Ha ( R k\ty < \\ u n N — u\\H a (B(0,N)\Cl) + 
\\ u n N \\H°>(R k \B(0,N)) + \\ u \\H a (R k \B(0,N)) ■ 

These last two terms converge to 0, and 

\\Un N — u\\H a (B(0,N)\n) < \\ u n N ~ ^A?"||if"(B(0,Ar)\Q) + \ \un ~ u\\H a (B(0,N)\n) 

X 

~ iV ~^ ~ M l|H Q (IR 3 \n)- 

This finishes the proof. □ 
Now, we can formulate the main result. 

Theorem 5. Let us assume u G D(A(0)) and g G C 1,v ((0, T); M, N+2 ), then there exists a 
solution of the problem 

^ = H{z(t))u+ H(z(t), z{t)) V« + (H(t) + H(t)V - |) G, 

B (t)u = 0, 

ii = ~Ki(z(t)) + R „ 2 } q ,., [ {u(x) + G)j{x,z(t))a"H 2 (x), 

u(0) = u o ,z(0) = 0. 

such that u G C a ([0,T); H s (tt c )) and z G (C 1,Q! ([0, T))) N+2 . Where I < s < \ + e and 
a < mm(l — (| + e) _1 s, i]). 
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For the sake of simplicity of notation we shall write L 2 = L 2 {Vt c ) and H s = H s (tt c 
Proof. We will apply the method of successive approximation. We set 

u°(t) = m , t e [0, oo) 
and let z® be a solution of the following problem 

(u°(x) + G)j(x,z(t))dH 2 (x), 



dt Zi ° pM z °^ + m 2 (s i (z))j Si 



(0) 



z? = 0, 



on the maximal interval of existence [0, To). Next, we define u n+1 as an unique solution (on 
interval [0, T n )) of the following problem 



du 



n+l 



H(t)u n+1 + H(t)Vu 



dt 



n+l 



Hit) + H{t) V -^\g 

B G u n+1 = 
u n+1 (0)=u Q . 



(12) 



Subsequently, we define z™ +l as a unique solution of 



d 

dt' 



1 



n+1 - x -^(z n+1 (t)) + 



l 



m 2 (s t ( z )) j Si 



[ (u n+ \x) +G)j(x,z n+ \t))dH 2 (x), (13) 

JSi(Q) 







on the maximal interval of existence [0,T n+ i), where T„ + i < T n . 

We shall show that sequences u n and ^ n are equibounded and equicontinuous. 

We can solve the parabolic equation ( fl2l using methods from previous section. Let as 
denote H(t) = H(z(t)) and H(t) = H(z(t), z{t)). 

We will show that the equation ( fl2l) has a unique solution. Indeed, let us denote by T an 
operator defined as follows 



G da. 



Tw[t) = U(t, 0)u + j U{t, a) (^5w{a) + H{a)Vw{a) + (h(o) + H{a)V - ^ 

We will show that the operator T possesses exactly one fixed point. First of all we shall show 
that there exist T n > and e n > such that 

B{u ,e n ) C C{[0,T n ];H s ) 
T : B(u ,e n ) -> B(u ,e n ). 



Let us notice 



sup \\H(s)\\ Loom < 

s6[0,T„_!] 



sup 

(2,t;)eB(eo)xS( 6 (,) 

+ sup 

se[o,Ti] 



(is) ( V ^l' a; 2' :C 3) A ^ 1 ' X2 ' X3 ' S )) 1 

/ (V( Xlja . 2 ^ 8 )A(ari,ar 2 ,a;3,s)) 1 



L 00 (R3 X R + ) 



+ 

<M l5 
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where B(e ) and B(e' Q ) are balls in C 1 and C a respectively, such that the solution z° is con- 
tained in B(e ) and if in B(e' ). 

Since G has a compact support it is not hard to see that 



sup 

se[o,T n _i] 



H(a)+H{a)V-^)G 



< M 2 . 



L 2 (f2 c ) 



Using the above inequalities and Corollary [4] we can write 



Tw(t)-u \\ Ha <c / L4(<r)C/(a,0)M 



c / ||Z7(t, tr)jy(o-)Vw(o-) || [z?(A(t))jia]/s do- + c 



[D(A(*)),L 2 ]p 
t 



da + c \\U(t, a)5w(a) I 



5 



U(t,a)[H(a) + H(a)V- — )G 



da. 



[D(A(t)),L 2 ]f3 



Using the methods from [23] one can show that the following inequalities 

C II II 



\A(t)U(t,a)u\ 



< 



u\ 



[D(A(t)),L 2 ]g ~ (t - a) X +P- 6 » »[J>{Mfr))M] e 



\\U{t,a)u\\ [D{A{t))Mh < \\u\\ L2 
hold for u G D(A(a)). Next, we obtain 

\\Tw(t)-u \\ HS <c^ -^\\ Uo \\ [D{m)M] da + c5^ 



\w{a)\\ L da + 



o {t-af 



\H(a)Vw(a)\\ L2 da + c 



H(a) + H(a)V 



d_ 



G 



da < 



l 2 



cT^\\u \\ [D(AmM]e + cT^ ( (5 + Mi) sup \\w(a)\\ HS + M 2 
where 6 > (3. But w G B(uo,e n ), hence, we choose T n and e n such that 

cT^ ((5 + M^wia)^ + M 2 ) + e n T n 1 -^(«J + M x c) < e n . 
Now, we show that the operator T is a contraction. Indeed, 



+ 



(14) 



\Tu(t)-Tv(t)\\ Ha <c5 / \\U(t,a)(u(a)-v(a)) 



\[D(A(t)),L 3 ] p 



da + 



c \\U(t,a)H(a)(X7u(a)-Vv(a)) 



[D(A(t)),L 2 h 



da < cS 



o {t-aY 



\u{a)-v{a)\\ L da + 



(t - *)P 



\H(a)(Vu(a) -Vv{a))\\ L2 da < (c5 + cM 1 )T n 1 " /3 sup \\u(a) - v(a) 

^e[o,T„] 



\H S ' 



We find e n and T n such that T is a contraction and the inequality ( TT4l ) is satisfied. Hence, 
from Banach fixed point theorem we obtain that T has a unique fixed point. Let us notice 
that e n = 62 and T n = T 2 for n > 2. 
We will show the following 
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Lemma 7. The sequence u n is equibounded and equicontinuous on [0, T 2 ]. 

Proof. Equiboundedness follows from previous considerations and construction of the sequence 
u n . Let us denote by a an upper bound of u n , i.e. 

sup ||ii n (cr)||H s < a. 

<tG[0,T 2 ] 

Now, we show that the sequence u n is equicontinuous. Namely, 



+ 



u n+ \t)-u n+ \r)\\ HS < \\U(t,0,)u o -U(r,0)u o \\ H s + 

d , ,. , 

da + 



(U(t, a) - U(t, a)) ( 5u n+1 (a) + H(a)Vu n+1 (a) + ( H{a) + H(a)V ~ t^) G 



U(t, a) [ 5u n+ \a) + H{a)Vu n+1 {a) + ( H(a) + H(a)V - i~ I G 



n+l 



+ 



It is not hard to show the following estimate 

h < c\\u \\[D(A(p)),L 2 ] e \t - r| e_/3 . 

Next, 



d 



da 



da 



h + I 2 + h- 



t r t 



OJt 



A(,i)V(,i.a) (r)V'+V) - H(a)Vu" +i -(a)+ ( H(a) + H{a)V - ^ ) G 



[D{A{r,))Mh 



t rt 



Cl 



5\\u n+ \a)\\ L , + \\H{a)Vu n+ \a)\\ L ,+ 



n+l. 



H(a) + H(a W - 



d_ 
da~ 



G 



h 
drjda 

drjda 



((5 + M l )a + M 2 )c 2 \t-T\ 1 



Subsequently, we estimate the last term J 3 : 



h < c 



Cl 



U(t.a) (rtV'+V) + f/((T)V//" +J -((T) - ( H(a) + H(a)V - ^- ) G 



da < 



(t - ay 



S\\u n+1 (a)\\ L 2 + \\H{a)Vu n+1 {a)\\ L , + || ( H(a) + H(a)V - ) G\\ L2 ) da < 



,n+l 



[D(A(t)),L 2 ] fJ/ 

d 



da 



((5 + Mi)a + M 2 ) c 3 \t - t 



1-/3 



Finally, 



\u n+ \t) - u n+1 {T)\\ H s <C(\t- r\ e -P + \t- r] 1 -' 3 ) . 



□ 



Now, we show the following 
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Lemma 8. The sequence u n has the following properties 

\\u n \\H°{m?\B{td,N)) — ► 
as N — > oo uniformly with respect to n. 

Proof. Let us denote by B C N = R 3 \ B(0, N). First of all we define the smooth map §n : Q 
[0, 1] as follows 

' 1 if r > iV 



ifr<f, 



$jv(r) 

where r 2 = x\ + x| + x\. We require that function <3>at fulfill the following condition 

l* (fe V)l<^. 

Subsequently, we define the map := ^^u n . One can easily show that u% solves the problem 

= Hu n N + HV($ N u n ) - u n HV$ N + $ n Ih + HV- — )G- 

-H<$> N u n - 2(H - f)V$ N (H - f)Vu n 

with initial condition u%(0) = Uq = &nUo- It is enough to show that ||Wjvll-ff s (n c ) — * as 
N — > oo uniformly with respect to n. 
Let us notice that 

\\U(t,0)uo\\ H *((ic) < \\uo\\h°(cic) + ciT® \\uo\\[D(A(o)),L 2 ] e (si°)- 
Now, we obtain the following estimate 

II n I Ml ^ ll iV|| i rri6—B\\ N\\ i 

SUp \\U N {a)\\ H s^c) < C\\U \\h°(CI c ) + CJ-2 \\ U \\[D(A(0)),L 2 \b{^ c ) + 

^e[o,T 2 ] 

+cT 2 1 " /3 ( sup |K(<7)||H. (n c)+ sup || ( H(a) + H(a) V - ) G\\v { & =) J + 

We[0,T 2 ] <re[0,T 2 ] V 0°~ J -r ) 



Hence, 



sup \\u%(a)\\ H ^n^ < c\\u ||//«(n c ) + c|jn \\[D{A{o)),L 2 ] e {n c ) + 
^e[o,r 2 ] 

+c sup || (//( CT ) + ^(a)V-|-)G'|| L2(B c )+ (-1 + 1) sup ||« n (<T)|| H . (B « =,)■ 
<re[o,T 2 ] V oa J f \iV z iV/ CTe[0iT2 ] T 

From above inequality and sup CTe [ 0T2 ] ||m"((t)||h s < a the lemma follows. □ 

Now, we examine the properties of the sequence z n . Let us recall that z n+l is a solution 
of the problem 



^<" +1 4 K,(2 " +1(t)) + ^^/ a( o,( G + M " +1W )^ 2 " +1(t))dW2W ' 

^ n+1 = 0. 
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First of all we show that the above problem has a unique solution. We define the operator T 
as follows 

Tz.it) = 1 jT KiizWds + jT fowls.^ I (G + u n +\x))j(x,z)(t))dH 2 (x)ds. 

We will show that there exist e n+ i and T n+ i such that 

B(0,e n+1 ) cC([0,T n+1 );M JV+2 ) 
.F : 5(0, e n+1 ) -> 5(0, 

and F is a contraction on B(0,e n+ i). Let us compute 



1 /" T " +1 X /" 

|^(t)-0|<- / K i (z( S ))ds+- / 

P« JO Pi JO 



ln+1 



A 7 n^iz)) J s 



(0) 



w n+1 (x) + G dH\x)ds 



< 



, V c 1 + T 7^(^> +-6), 



where 



Ci = SUp SUp |/Cj(z)|, 
l<j<Af+2ze_B(0,£ 2 ) 

/3 = inf pi, 

l<i<N+2 

K— sup sup \J(x,z)\, 

x£dn c z£B(0,e 2 ) 

k= inf inf H 2 (SAzj), 

l<i<N+2 zeB(0,e 2 ) 



b = sup sup / \G\dH (x) 

l<i<N+2 zeB(0,e 2 ) JSi(0) 

and 7 is the constant from the trace theorem. 

Now, we show that T is contraction on B(0, £ n +2)- 

1 f^n+l 



Tz(t) - T~z{t) < - / Ki(z(s)) - iu(z(s)) 
P Jo 



+ 



r T n+i r 
Jo JSi(O) 



u n+ \x) + G) 



J(z(s)) J(z(s)) 



n 2 {Si{z)) H 2 {Si{z)) 



ds + 



< 



\{c + kc 2 r ya + c 2 b)T n+1 sup \z(s) - z(s)\. 
P se[o,T n+1 ] 



Hence, form contraction principle the solvability of our problem follows. Let us notice that 
e n = 62 and T n = T 2 for n > 2. 
We have the following 

Lemma 9. The sequences z n and z n are equibounded and equicontinuous. 
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Proof. From above consideration follows that z n is equibounded. We show that it is equicon- 
tinuous. Namely, 



z n+1 (t)-z n+i {a) 



n+1, 



< 



1 





Ki(z(s)) 



ds + 



pj a n\s i {z^))j Si 



(0) 



u n+1 (x)J(x,z)(t) dH\x)ds + 



+- 



(0) 



GJ(x, z){t) dH\x)ds < - ci + —ia + — b )\t - a 



Subsequently, we show that z n+ % is equibounded, 



z n+l \ < i [ct + '—^a+'—b] . 



1 







K 







K. 



K 



K 



k 



k 



k 



Next, equicontinuouity follows from the estimate 



z n+1 (t) - z n+ \a)\ < - \^{z n+ \t)) - Ki (z n+1 (a))\ + 



1 

Jsifp) 







dH 2 < c\t - o-l 



n 2 (Si(z^)) v v " H 2 (Si(z^)) 

where we applied Lipschitz continuity of z and Holder continuity of u n+1 . 

From the above three lemmas, Lemma [6] and Ascoli Theorem the proof follows. 



□ 
□ 



6 Summary 

In this paper we have considered the system modelling evolution of crystals 3-D from supersat- 
urated vapor. We have shown local in time existence of solutions. Let us mention that we could 
not expect global in time existence of solutions (topological catastrophe could occur). Let us 
stress that uniqueness of solutions is still an open problem. In the forthcoming article (see 
[T5] ) we will show local in time existence of solutions to the problem with Stefan conditions. 
In this problem we shall apply Galerkin method and we shall obtain less regular solutions. 
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7 Appendix 

This section contains an auxiliary and technical result. 
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Lemma 10. Let us assume that Q is an admissible shape, i.e. Q = W x [—L,L]. Then there 



exist continuous maps hi (i = 1,2, N + 2) with compact support such that an 

exist a.e. and §|, e L 2 {Vt c ) (i = 1,2, JV + 2 and fc, / = 1,2,3/ Moreover the following 

is true 

dh. 



5,(0) 



= 5. 



Proof. It is enough to construct maps hi for i — 2,T, B. 

Let us introduce TV + 2 cut-off maps e C^°(IR 3 ), < rji < 1. The following positions of 
lateral facets are possible (see Figure 7) 



V 2 




1 


*x 2 


V 2 | 






S 2 \ x 




sA 

-\ — 




X 





Figure 7 

Let us denote 5 = 10diam(Q). For i — 1, 2, N we define cut-off function as follows 

if (xi,x 2 ) ^ (Li) 4e or |(xi,x 2 ,x 3 )| > 5 



rii(x 1 ,x 2 ,x 3 ) 



1 if (xi,x 2 ) G (Lj) e and |x 3 | < f, 



^(^1,2:2,2:3) 



where e is sufficiently small. Subsequently 

if x 3 > or \(x\, x 2 , x 3 )\ > 5 

1 if x 3 G [L - £, L + e] and | (xi, x 2 ) \ < f 

and t)b(xi, x 2 , x 3 ) = t)t(xi, x 2 , — x 3 ). Next we denote 

fyr(xi,X 2 ,X 3 ) = 77t(^1,^2,^3)^3, 

h B (xi,x 2 , x 3 ) = i] B (x 1 , x 2 , x 3 )x 3 . 

and 

{x 2 r] 2 (x 1 ,x 2 , x 3 ) if (xi, x 2 ) E A 2 

x 2 rj 2 (x!,x 2 , x 3 )0 2 (xi) if (x u x 2 ) e B 2 
x 2 7] 2 (xi,x 2 , x 3 )(/) 2 (xi) if (xi, x 2 ) G C 2 , 
where the sets A 2 ,£? 2 and C 2 are explained in Figure 8. 




4e 




Figure 8 
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where 02 (#i) and 2 (^i) are chosen in such a way that the required properties are fulfilled. 

□ 
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